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During the grant period starting July 1, 1994, our major effort has been on the following 
two problems: 

1. Temporal behavior of heavy Oxygen ion outflow in response to a 
transverse heating event. 

2. Continued effort on ion heating by lower hybrid waves. 

We briefly describe here the research performed under these topics. 

1. Temporal response of 0 + ions to a transverse heating event 

Observations of superthermal 0 + ions from Dynamic Explorer 1 have shown both up- 
flowing and down-flowing populations A more interesting observation from AKEBONO satellite 
is that a cold ( < 1 eV ) population of 0 + ions is commonly found at altitudes of several thousands 

of kilometers. Knowing that the 0 + ions originate in the ionosphere where they are generally 
cold (< 1 eV) and gravitationally bound, the occurrence of cold 0 + ions at high altitudes is a 

mystery. Our recent modeling efforts [Singh, 1 995] explain that the creation of the superthermal 
up- and down-flowing ions along with the generation of cold 0 + ions at high altitudes are all 
related phenomenon Using a mesoscale particle-in-cell code, we investigated the time-dependent 
flow of 0 + ions in response to a brief (a few minutes) transverse heating event at altitudes of 
about 2000 km. The heating was chosen to produce ions up to about 15 eV in perpendicular 



energy. As expected, the ions having energies > 8 eV escape, while those with lesser energies are 
gravitationally bound, but they rise to high altitudes depending on their initial energy. Such ions 
become progressively colder as they move upward because their kinetic energy is converted into 
gravitational potential energy Near the apex of their trajectory, the ions are quite cold. During 
their return journey downward, the ions are again accelerated to relatively large downward 
velocities as observed from the Dynamic Explorer. 

The falling ions mirror upward when they reach the ionospheric altitudes. Thus the 
initially heated ions are trapped inside the flux tube and undergo a repeated bounce motion. Since 
initially ions have a wide range of perpendicular energies after the heating, the ions cover a wide 
range of altitude as they execute their bounce motions The result of this bounce motion is that 
the ions become progressively colder with increasing altitude with the coldest ions occurring at 
the highest altitudes considered in the model 

The rising and falling superthermal ions and the very cold 0 + ions at high altitudes seen in 
the model have characteristics similar to those reported from observations from the Dynamic 
Explorer and AKEBONO satellites. The details of the models and results can be found in the 
publications enclosed along with this report. 

2. Ion heating by lower hybrid waves 

Lower hybrid waves are known to heat ions. In the auroral plasma lower hybrid waves 
are a part of the broadband VLF waves driven by auroral electron beams carrying both upward 
and downward currents. Normally the waves driven by the electron beams have wavelengths so 
long that they do not affect the ions. In other words, this implies that ions do not resonate 
directly with the fast waves driven by the electron beam. It has been suggested that lower hybrid 
waves evolve nonlinearly into short wavelength slow waves When the waves have sufficiently 
evolved, they might resonate with the ions and heat them perpendicularly. 

During the last six months we have attempted to simulate this nonlinear evolution on a 
computer. The computer model uses a 3-D PIC code, as demanded by the nature of the problem. 



We have found that a sequential code cannot serve the purpose This has led to develop a parallel 
code. The development of the parallel code is being continued now under a new grant (NAGW- 
4470). 
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ionosphere, J. Geophys. Res., in press, 1995 

2. N. Singh, Fountain-like flow of heavy oxygen ions from the earth's ionosphere in response to 
transverse heating, IEEE Trans. Plasma Science, Special Issue on Images in Plasma Science, 
reviewed and resubmitted, Nov. 1995. 



Time Response of 0 + to a Weak 
Transverse Ion Heating Event in the Polar Ionosphere 


Nagendra Singh 


Department of Electrical and Computer Engineering 
University of Alabama 
Huntsville, AL 35899 


Submitted to JGR, May 17, 1995 



Abstract: 

Time response of 0 + ions to transverse ion heating in the polar ionosphere is studied by 
means of a particle code. Here we deal with relatively low levels of heating over a few minutes 
yielding superthermal ions energized up to a few eV. Under the influence of the upward mirror 
and downward gravitational forces, most of the heated ions are trapped between the topside 
ionosphere and some high altitude depending on their perpendicular energies. A few ions with 
sufficiently large perpendicular energy escape the gravitational pull, and the flux of such ions 
increases with the heating level. The trapped ions bounce back and forth. During the first 
upward transit of the heated ions, the transient outflux from the topside ionosphere is as large as 
10 ions cm / s and it subsequently decays to quite small values in a few tens of minutes. A 
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large transient outflux (-3x10 ions cm ~ / s ) consisting of unheated ions develops even 
below the heating region lasting over just a few minutes. As the ions having relatively large 
energies continue their upward joumev, and the less energetic ones begin to fall down, the 0 + 
drift velocity ranges from a relatively large negative value (downward) just above the dense 0 + 
plasma in the topside ionosphere to a large positive value at high altitudes. Substantial flux of 
downward moving ions is also seen. During subsequent upward motions, the bouncing ions form 
a well-defined, quite sharp expansion front. The bounce motion generates a succession of plateau 
formation over a period of several hours in response to a brief heating period of just a few 
minutes. This leads to nearly periodic oscillations in local densities, with the oscillation period 
increasing with altitude. In a time of about one day, the oscillations in the density cease due to 
the spatial dispersion of the heated ions having differing initial energies. Eventually the 0 + 
density profile becomes nearly stable, but remains extended with significantly enhanced 0 + 
density at high altitudes. A consequence of trapping of the ions by the downward gravitational 
force is that the kinetic energy of the ions decreases with the increasing altitude; the coldest ions 
appear near the top of apex of the trapped ion trajectory while the wannest ions appear near the 
bottom of the flux tube. Such cold ions are likely to contribute to the thermal 0 + population 
observed at altitudes of several thousand kilometers. 
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Introduction: 

Transverse ion heating is a common feature of the Earth's magnetosphere. Such heating 
occurs over a wide range of altitude and latitude, and it is commonly invoked to explain the 
escape of gravitationally bound heavy ionospheric ions such as 0 + [Lockwood et al, 1985; Li et 
al, 1988; Waite et al, 1985], In order to escape the Earth's gravity, 0 + ions in the topside 
ionosphere should be energized to about 10 eV. Escaping 0 + ions with energy in the range from 
a few hundreds of eV to several keV are commonly found at high altitudes in the auroral plasma. 
Such ions flow upward with conical pitch angle distributions, which are created when 
transversely heated ions are acted upon by the mirror force, and the cone angle progressively 
decreases with altitude. Such considerations for the measured energetic 0 + ions show that the 
heating must occur at a high altitude above about 5,000 km, where abundance of 0 + ions is 
normally too small to generate the observed large flux of about 10 ions cm s . Therefore, it 

is logical to ask how the gravitationally bound 0 + ions are transported to such high altitudes. 

Barakat and Schunk [1983] suggested that enhanced electron temperature ( 1 0 4 £ ) can 
facilitate the escape of 0 + ions; the escape occurs by enhanced polarization electric fields. Singh 
et al [1989] showed that the enhanced electric field associated with the density gradients in the 
auroral cavity can lead to a relatively large 0 + flux. Li et al [1988] showed enhanced 0 + flux can 
be generated by raising parallel and/or perpendicular temperatures of 0 + ions at some low 
altitudes. Using isotropic hydrodynamic model, Gombosi and Killeen [1984] and Schunk and 
Sojka [1989] showed that the low altitude frictional heating leads to the upflow of 0 + ions. The 
frictional heating of 0 + ions and associated dynamics were treated by Korosmezey et al [1992] 
using a more sophisticated hydrodynamic model, in which temperature anisotropy was allowed. 
Recently Wilson [1994] pursued a similar idea using a semikinetic model. Lorcmc and St. 
Maurice [1994] also studied the temporal behavior of 0 + outflow in response to a fractional 
heating at altitudes below 500 km using a semianalyticai kinetic model in which the polarization 
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electric field was assumed to remain uniform and constant. Since the frictional heating heats 0 + 
to about ~! eV, the heated ions do not escape; after an initial uptlow the ions fall under the 
influence of the gravity [Lorcmc and St. Maurice , 1995], 

The purpose of this study is to examine the dynamics of the non-escaping ions, which are 
energized to a few eV and remain gravitationally bound. Such ions show an interesting set of 
dynamical features, which have not been studied before. These features include (1) simultaneous 
occurrence of down-flowing 0 + below certain altitude, and up-flowing 0 + above it during the 
early phase of the evolution; (2) the down-flowing 0 + eventually bounce back and forth as they 
are trapped between the bottom of the flux tube and some high altitude depending on their initial 
kinetic energy; (3) when the trapped ions move upward, they create a plateau in the density; the 
leading edge of the plateau is like an expansion front, which moves with a nearly constant 
velocity; (4) recurring formation of the plateau is seen; this creates oscillation in local densities 
with a periodicity which increases with increasing altitudes from about 2 hours at a low altitude 
(geocentric distance r = \0 4 km) to several hours at high altitudes, and (5) eventually the 
oscillations subside as the ions with different energy gains during the heating period disperse 
giving rise to an 0 + density distribution substantially extended in altitude compared to the one 
without the heating. The top of the extended distribution is near the apex of the trapped ion 
trajectory for the ion with the maximum initial energy. Due to the gravity, the kinetic energy of 
the trapped ions decreases with altitude producing very cold ions near the top of the density 
distribution. 

Recently Chandler [1995] reported observations of the down flowing 0 + . The downward 
velocity of the ions was typically <1-5 km/ s below the altitude of 4000 km . Our simulation 

shows how a brief transverse ion heating leads to the creation of a population of falling ions 
during the very early phase of the bounce motion. The evolution of the falling ions from the 
initially upflowing population of heated ions is critically examined. In connection with the cleft 
ion fountain [Lockwood and Horwitz , 1985], Horwitz [1984] and Horwitz and Lockwood [1985] 
studied 2-dimensional trajectories of bouncing 0 + superthermal ions. They found a "hopping" 
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trajectory in r - X plane where r is the geocentric distance and X is the latitudinal angle. The 
conversion of bounce motion into a hopping trajectory is caused by the ExB drift where E is 

the convection electric field. In the present paper, we have not included the convection field, as 
our main goal is to understand the basic dynamical features driven by the brief period of 
transverse ion heating. We report here the time constants of the different phases of the evolution 
from ion heating to the generation of extended 0 + density profiles. We also report here both the 
downflux and upward flux of 0 + ions generated by the transient heating and compare the latter 
with those produced in the steady-state models. The temporal and spatial evolutions of the fluxes 
are discussed. 

Using data from the Akebono satellite, [Abe et al, 1993; Peterson et al 1993] reported a 
common occurrence of a stationary, dense cold 0 + population extending to altitudes of several 
thousand kilometers. How such cold 0 + ions are transported to such high altitudes has remained 
a mystery. We demonstrate here that the gravity in combination with transverse ion heating at 
low altitudes producing superthermal 0 + ions creates a cold 0 + population extending to high 
altitudes as observed 

The rest of the paper is organized as follows. In section 2 we briefly describe the model. 
In section 3, we discuss the dynamics of the superthermaily heated ions by the process of 
transverse ion heating. 0 + fluxes associated with the up- and down-flowing ions are discussed in 
section 4. Conclusions and discussion are given in section 5. 

2. Numerical Model 

We study plasma flow in an open flux tube as shown in Fig. 1. The topside ionospheric 
base is assumed to be at r = r b . We simulate the flux tube up to a geocentric distance r ~ r max • 
We assume that Earth's magnetic field B(r ) s B 0 {r b )(r b / r) 3 , where B a (r b ) is the magnetic field 
at r -r b . In order to study the plasma flow, we employ a particle code [ Wilson et al , 1990], The 
model treats a collisionless plasma. Therefore, the boundary at r = r b is taken to be at a 
sufficiently high geocentric distance r b = 7350 km or altitude // = 980 km. Note that in this 
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paper (height) /j = r-Re, where Re is the Earth's radius. We inject macroparticles 
corresponding to 0 + and H + ions into the flux tube at the boundary at r = r b and subsequently 
follow their dynamics under the influence of the gravity, polarization electric field and the mirror 
force associated with the gradient in the geomagnetic field. The injected particles are chosen from 
Maxwellian velocity distributions with a zero drift for 0 + and a non-zero drift for H + The choice 
of H + drift does not significantly affect the outflow of 0 + ions in response to the heating. In this 
paper we present results for a drift velocity of 1.55 Km, where V tH is the thermal velocity 
corresponding to the temperature T 0 of H + ions at the base. Only the charged particles with 
parallel velocity V\\>0 are injected into the flux tube. We assume that 0 + and electron 
temperatures at r = r b are also T 0 . When a steady-state polar wind flow is set up in the flux 
tube, we initiate the transverse heating event for 0 + ions. The heating is assumed to be extended 
and it occurs for all altitudes h > 1 940 km (r > 83 1 1 km) . The heating at such altitudes is 
expected to be caused by wave-particle interactions in contrast to the frictional heating occurring 
at low altitudes discussed earlier. The heating occurs in a random fashion as described by Brown 
et al [1991], Whether or not a particle is heated during a time step A t depends on the 
probability density given by 

P(A w i )= — ! — exp [-UA w, /a ) 2 ] (1) 

2 no 

where A h ) 1 is the impulse of energy given to a 0 + macroparticle in the degree of freedom 
perpendicular to the magnetic field. In implementing this heating scheme, a random number 
is chosen from a uniform distribution from 0 to 1 for each 0 + macroparticle lying in the desired 
altitude range. This number is compared against the cumulative probability associated with the 
above probability density and when the latter exceeds R„, the energy of the particle is 
incremented by A w L . Brown et al [1991] showed that the above scheme gives a heating rate 
given by 

= a / 1.14 A t ( 2 ) 

dt 
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where At is the time step used in advancing the motion of the 0 + macroparticles, k is the 
Boltzmann constant, and 7]_ is the average perpenaicuiar temperature of the heated ions. 
Equation (2) defines cr in terms of the heating rate. 

In our simulations we assumed that 0 + and H + densities at the base r = r b are 
= 0.9 n 0 and n^ + = 0.1 n Q , where n Q is the total plasma density there. At the base H + 
flow was assumed to be supersonic with an average flow velocity u (r = r b ) = 1.55 V tHi where 
V sH = kT 0 Tntjj + with m H+ being the mass of an H + . For 0 + we assumed u Q+ (r = r b ) = 0. 
Time step in advancing the motion of the macroparticles was A t = 1 1 s . The polarization electric 
fields were calculated by assuming that electrons obey the Boltzmann law with a constant 
temperature T a and the plasma remains quasineutral; this yields eiectric field 
E = -(kT a / e)xn~ l d n e / d x , where n e = n Q . • In the following discussion we use the 
following definitions and normalizations: density n=n/n Q , velocity V - V IV tH , electric field 
E -E / E a , and flux r = r l r o where n and V denote the density and velocity of charged 
particles, respectively; r o = n Q and E 0 - k B T Q l eL, with L - V^j / g Q , which is about 
3000 km for T o =0AeV and g a = 9.8 ms~ 2 and E 0 = \jjVlm. For the chosen parameters 
of our calculations ( n Q = 10 4 cm~ 3 and T a =0.3 eV) V tH =5.5 km/ s, E a = 1 uV / m and 
r n - 5.5 x 10 9 ions cmT 2 s ~ 1 . 

3. Numerical Results 

In this section we first briefly describe the ambient polar wind, followed by a detailed 
discussion on the transient outflow of 0 + ions, including interpretation of the numerical results. 

3.1 Ambient Polar Wind 

Starting with an empty flux tube, a polar wind consisting of 0 + and H + ions is set up by 
injecting charged particles into the flux tube at r = r b (Figure 1). It takes several hours to 
establish a stable polar wind as determined by the transit time of the slower 0 + ions, as they move 
upward and then fall under the influence of the gravity. Our numerical calculations showed that a 


6 



stable polar wind is set up in about 6 hours, which is rougniy the time for 0 + ions with a velocity 
of ~ 1 kmi s to travel from the boundary at r - r h to an altitude r - 10 4 km and then fall back 
to r = r b . Figures 2a to 2f show the state of the ambient oolar wind. Since the state shown in 
these figures serves as initial condition for further calculations on 0 + dynamics in response to 
transverse heating, we label this state by t = 0 . Figures 2a and 2b show the phase space plots for 
0 + ions in r-V \ and r-V L planes, respectively; v\ \ ana \\ are the components of velocity of 
charged particles parallel and perpendicular to the magnetic field. The corresponding plots for H + 
ions are in Figures 2c and 2d. 0 + ions are almost in static equilibrium, except for a few ions 
having initially large perpendicular velocities. These ions are oropelled to relatively high altitudes 
under the influence of the mirror force. However, with the temperature T a - 0.3 eV even these 
ions do not escape and are eventually turned downward. The 0 + phase space in r-V\ \ is quite 
symmetric about V\\ = 0 indicating an average drift velocity of zero. 

Figure 2c shows that the light H + ions flow out and attain relatively large supersonic flow 
velocities all along the flux tube. Figure 2d shows the expected cooling in the perpendicular 
degree of freedom as H + ions move upward and eventually escape from the top of the flux tube. 

The density profiles of the 0 + and H + ions in the ambient polar wind are shown in Figure 
2e. The crossover point occurs at r = 9\00 km . ~his figure also shows the normalized 
polarization electric field E associated with the ambient Doiar wind. The maximum electric field 
occurs at the base (r = r b ) and its volume is E = 9.6 E 0 £ 0.96 Im . The electric field profile 
is determined by the two density profiles in the figure because E oc n~ x d n e t d x where 

"*=v + v- 

Profiles of the average flow velocities of 0 + ana H + ions are shown in Figure 2f. For 0 + 
the average drift velocity is nearly zero all along the flux tube, as expected for a static equilibrium. 
On the other hand, H + ions are nearly supersonic all along the flux tube; near the base the average 
drift velocity is 1.55 . There is no net outflux of O' - ions in the ambient polar wind. The 

outflux of H + ions is shown in the figure; we have plotted here F which is the normalized 
outflux defined by F (r) = r / (r) u ;r (r) (r /r h )\ First of all we note that F ( r ) 
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is nearly constant as a function of r, barring some fluctuation inherent to a PIC code. The 
constancy of r indicates that the code is fairly accurate. The constant value of r is 
determined by the boundary condition and it is r = 0.155, which amounts to a flux 
r =0.155 n 0 Vtfj s8.5x 10 ions cm ~s at r = r b . This flux is determined by the boundary 
conditions on H + average velocity (1.55 V dl ) and density (« /r = 0.1 n a ) at r = r b . 

3.2 Time Response to Transverse Heating 

Our main goal here is to study the temporal behavior of transversely heated ions at 
relatively low altitudes. We consider here only a low level of heating which produces 
superthermai ions in a few eV energy range. We examine the temporal behavior of (a) the phase 
space in r -V\\ and r-V L planes, (fi) 0 + density and average flow velocity profiles and (c) 
temporal variation of the local ion densities. Figures 3 a and 3b show the phase space plots after 5 
minutes of heating with a heating rate d Id t (kT^) = 0.04 kT 0 / s , which gives cr = 0.5 eF from 
equation (2). At this time we note that the perpendicular velocity distribution of the ions is 
significantly modified by the perpendicular acceleration of ions up to V L = 2.5 V t „ and maximum 
perpendicular energy W, max =51 kT 0 . Assuming kT Q = 0.3 eV , = 15 eV . 

For the flux tube we consider, any ion exiting from the top of the flux tube is essentially 
lost or escaped for our purposes. Ignoring parallel electric field, we can estimate the total energy 
required for such an escape and it is 

r=r max i | 

»'»«- J « m(r g 0 RU ) (3) 

r=?j, r b Tnax 

where g Q is the gravity at r = Re , the Earth's radius. 

For the parameters of our calculation (r b = 7350 km, /- max = 31350 Awi, g 0 =9!&ms~ 2 . 

Re = o371 km ) , lV es = 5.7 eV or 19 kT 0 . Therefore, we expect that heated ions with 
1/2 

V ]_ > (2 W cs / m Q +y = 1.5 V , H will eventually escape from the flux tube. We point out that if the 
top boundary were taken to be at infinity, i.e., — >• oo, W es = 7.5 eV and perpendicular escape 

velocity is then only slightly larger. 
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top boundary were taken to be at infinity, i.e., r raax — » oo, IV CS = 7.5 eV and perpendicular escape 
velocity is then only slightly larger. 

From the above consideration, we expect that some of the heated ions in Figure 3b will 
eventually escape from the flux tube as the mirror force accelerates them upward. On the other 
hand, those ions with V ± < 1.5 V lH will be eventually returned downward after reaching their 
apex. This is clearly brought out from the phase space plots at later times as shown in Figures 3c 
and 3d for f = 20min., Figures 3e and 3f for / = 60min., and Figures 3g and 3h for 
/ = 120 min. At / = 20 min. , ions show parallel acceleration by the mirror force. When 
t = 1 hrs., the accelerated ions have reached the top of the simulated flux tube and their parallel 
average flow velocity is U\\ = 1.36 V tH - 7.4 kmts. Upflowing superthermal 0 + ions with this 
magnitude of parallel flow velocity have been observed from Dynamic Explorer 1 at geocentric 
attitudes shown in Figure 3 [Waite et al, 1985]. The phase space plots in Figure 3e shows that 
the average flow velocity of the upflowing ions varies almost linearly with the altitude. As the 
time passes, the heated ions having relatively low perpendicular energies begin to move down 
first. Figure 3g shows that at t = 120 min. ions are already falling (Fj ( < 0) at geocentric altitudes 
below r < 18000 km, and the falling ions are beam-like until they begin to penetrate the dense 0 + 
population at low altitudes. Thus, for heights h < 3000 km ( r < 10,000 km) it is difficult to 
distinguish falling ions from the ambient 0 + population, but above this height falling 0 + create a 
distinct beam-like population and maximum downward velocity in the height range of about 4000 
km is |F[|| < 0.5 V tH =2.5 kmis. At higher altitudes, the average downward flow velocity 
decreases with increasing height at t = 120 mins. At low altitudes below 4000 km also, the 
average flow velocity decreases and becomes quite smail, as the falling ions penetrate into the 
dense 0 + population. This yields a maximum flow velocity somewhere in the midaltitude range. 
Recently Chandler [1995] reported data on downflowing 0 + form DE-1 satellite; his data are 
confined to altitudes <4000 km in the polar topside ionosphere. The data show that the 
downward flow velocity of 0 + nearly monotonically increases with altitude and reaches a value 
~2 km! s at about 4000 km altitude. No observation is available on down flowing 0 + for 
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h > 4000 km . Our simulation shows that average flow velocity of 0 + ions should monotonically 
decrease with distance above a certain altitude, where it maximizes, but the height of the 
maximum downward drift depends on the phase of the evolution. We further discuss this in a 
later section on 0 + flux. 

The down-flowing 0 + event discussed above lasts for much longer time than that covered 
in Figure 3. Figure 4 shows the further evolution of the transversely heated ions including those 
flowing down. The downward flowing ions with flow velocity varying almost linearly with the 
altitude are seen in Figures 4a and 4c, which show the phase space plots in r -Vj | plane for 4, 

and 6 hours, respectively. Such downflowing ions recede from the top and eventually disappear 
as seen from Figure 4e for t = 8 hrs. However, they reappear as shown for t - 12 hrs. when 
the down flowing ions mirror and move upward to fall again, as discussed below. 

Beside falling 0 + ions. Figures 4a to 4h show another interesting feature of the 0 + 
dynamics in response to the brief heating lasting over 5 minutes. Figures 4a and 4b reveal that 
besides the bulk of 0 + below r = 10 4 km, a band of new 0 + cloud is appearing at higher 
altitudes as indicated by an arrow in Figure 4a. This band of 0 + ions in r - V\\ plane has a 
distinctive signature in the r-Fj_ plane as well, as indicated by an arrow. Near the cloud the 
density in the phase space is relatively high as marked by the dense packing of the dots 
representing ions. Figures 4c and 4d for t = 6 hrs. show that the 0 + cloud in the band move 
upward, and a new band appears to separate from the dense 0 + plasma at the bottom. Figures 4e 
and 4f for / - 8 hrs. show that both bands of 0 + cloud appearing in Figures 4c and 4d have 
moved up. At later times the ions in the first band begin to fall after reaching their apex, while the 
ions in the second band continue to move up; these features are seen from Figures 4g and 4h for 
/ = 12 hrs. These figures show the development of additional bands. The band structure 
appearing in the phase space plots of Figures 4a to 4h is a consequence of downflowing ions 
mirroring back and moving upward again. However, a question arises as to the cause for the 
more than one band to appear. The upward motion of the superthermai ions is self-modulating; as 
the ions move up modifying the density profile, the electric field is modulated in time. When the 
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in the electric field has an appreciable affect on the 0 + dynamics. This creates bunching of the 
superthermai ions; each bunch appears as a band in the r - V n phase space. We wiil further 

discuss the oscillation in the electric field after presenting the dynamic behavior of 0 + density 
profile 

Evolution of 0 + Density Profile: The repeated upward motion of the mirroring ions is 
more clearly seen by looking at the density and the average flow velocity profiles of 0 + . Figure 5 
shows the evolution of the 0 + density profile for times before t = 3 hrs., about the time when the 
mirroring ions begin to move upward for the first time after their first fall. For the purpose of 
comparison, we also give the 0 + density profile at t = 0 in the ambient polar wind. In response to 
the heating, 0 + ions move upward as they are accelerated upward by the mirror force, without 
any significant change in the flow of H + ions. By the time t = 1 hrs. the heated ions begin to 
reach the top of the simulated flux tube. As the heated upweiling ions disperse in the flux tube to 
higher and higher altitudes and some of them flow out from the top, the density beings to drop in 
the newly populated altitude region as seen from the plots for t > 90 min. The flux of 0 + in such 
an upweiling event will be discussed later. We notice from Figure 5 that at t = 1 80 min. , the 
density profile beings to expand again as evidenced by the formation of a ledge in the density 
profile below r = 14,000 km . This re-expansion is shown in Figure 6, which shows density 
profiles from / = 3 hrs. to 6 hrs. The re-expansion is accompanied by a well-defined expansion 
front as marked by arrows. The velocity of expansion front is nearly constant and it is about 600 
m/s. Locations of the expansion fronts correspond to the previously discussed formation of the 
first band of 0 + cloud seen in the phase space plots in Figures 4a and 4c. The evolution of the 
average flow velocity of 0 + for the expanding piasma in Figures-6 is shown in Figure 7. Note that 
the velocities shown are normalized to V* = 5.5 km! s . The location of the expansion front is 
where the velocity has a sharp jump from a small positive value to a relatively large negative 
value. The negative velocity above the locations of sharp jump corresponds to the falling ions, 
discussed earlier in connection with Figures 3 and 4. Thus, the expansion front separates the 
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upflowing ions just below it from the dowmlowing ions above it. Below the expansion front the 
average velocity is generally quite smail and it approaches zero in the dense 0 + plasma at low 
altitudes. However, at times the dowmlowing ions may dominate even below the expansion front 
as revealed by the average rlow velocity at t = 6 hrs. in Figure 7; at this time the average upward 
velocity near the expansion front is about 500 mis , while just below it the average downward 
velocity is about 1 km / s . Figure 6 shows that at t = 6 hrs. there is yet another new expansion 
front separating from the dense plasma at low altitudes; the new expansion front corresponds to 
the new band appearing in the phase space in Figure 4c. Further expansion of this new front is 
shown in Figure 8 by the short arrows. We notice both the old and new expansion fronts in this 
figure; the new front is at lower altitudes, while the old one has traveled to r = 2300 km by 
t = 8 hrs, as indicated by the long arrows. 

Electric Field Oscillations: The above dynamical behavior of the 0 + density modulates the 
electric field in time, which in turn affects the dynamics of the 0 + ions as they bounce back and 
forth. The inset of Figure 5 shows an example of oscillation in the electric field at r = 9230 km , 
a location close to the crossover point between the 0 + and H + density profiles in the ambient 
polar wind. Before the heating, the normalized electric field at this location has a constant value 
of about 3.2. In response to the heating, the electric field is seen to undergo oscillations, which 
are sensitively controlled by the slope of the 0 + density profile because in our model 
E = -(kT 0 !e)n; x dn e Idx, where n e ~ n Q+ + n JJJr . From Figure 1 we noted that the 
crossover point between 0 + and H + density profiles is at r =9100 km . The heating occurs for 
r > 83 1 1 km . Near and below the crossover point, the 0 + density is significant compared to the 
H* density. Therefore any modulation in the 0 + density profile is reflected in the modulation of 
the electric field. Far above the crossover point, where n Q+ « n } ^ , the modulation in the 0 + 
density profile does not by itself cause a modulation in the electric field. However, the 
modulation in the electric field at relatively low altitudes does perturb the H^ density resulting in 
some modulation in the electric field, even at high altitudes. 
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The inset of Figure 5 shows that in response to the heating, the electric field first sharply 
increases from its normalized value of 3.2 to 3.4 followed by a sharp drop and subsequent 
oscillations about a reduced levei of time average electric field. These temporal features are 
related to the changes in the slope of the 0 + density profile. However, due to the log scale used 
in the density plots in Figure 5, it is difficult to clearly see the changes in the slope near 
r =9230 km. When plotted on a linear scale in a localized region, the density profiles do show 
the change in slopes as borne out from the oscillations in the electric field. The oscillations in the 
electric field modulate the outflow of the 0 + ions forming the bunches seen in Figure 4. Even a 
relatively smail change in the electric field becomes significant in 0 + motion near the mirror 
points where the parallel velocity vanishes. In contrast, the FT 1- flow is only slightly affected by 
the oscillations in the electric field because the flow is highly supersonic for most parts of the flux 
tube. 

Oscillations in 0 + Density: The reformation of density expansion front in connection with 
the bounce motion of 0 + ions is found to cause well-defined periodic oscillations in the local 
densities. Figure 9 shows these oscillations in densities at some selected geocentric distances. 
The density oscillations at r = 9270 km , where the 0 + density is large, is relatively weak. 

However at higher altitudes the relative amplitudes of the oscillations become quite appreciable. 
The first peak in the oscillations corresponds to the first wave of heated ions as they move upward 
under the influence of the mirror force. On the other hand, the subsequent peaks are the 
consequences of the passage of the expansion fronts. The sharpness of the increase in the 
densities is explained by the sharpness of the expansion fronts. The density at any altitude slowly 
decreases after the passage of the expansion front because as the plasma cloud expands, it 
occupies a larger volume giving rise to a progressively decreasing density. Note that the 
oscillation periods progressively increase with the increasing altitude. This simply implies that it 
takes longer for an expansion front to reach higher altitudes. Note that the oscillations occur up 
to distances of a few tens of thousands of kilometers. 
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takes longer for an expansion front to reach higher altitudes. Note that the oscillations occur up 
to distances of a few tens of thousands of kilometers. 

The recurring reformation of the expansion front continues for a fairly long time until the 
superthermal ions disperse in altitude due to their different initial energies acquired during the 
heating period. We find that this time is about a day. The inset of Figure 10 shows the phase 
space plots of ions at / = 24 hrs. The ions are dispersed in altitudes up to r = 29000 km , where 

the apex of the trapped superthermal ions occur. The trapped ions are in a nearly diffusive 
equilibrium with average flow velocity approximately zero all along the flux tube. The density 
distribution associated with the phase space plots is also shown in Figure 10. For the purpose of 
comparison, we have also plotted the initial 0 + density profile just before the heating. The result 
of the heating is that the superthermal ions trapped in the flux tube significantly enhance the 0 + 
density up to altitudes of several thousands of kilometers. Such an enhancement may be part of 
the solution for the supply of 0 + ions to high altitudes where transverse ion heating produces 
energetic ions in keV energy range. 

The phase space plots in Figure 10 reveals an interesting feature of the trapped 0 + . The 
largest kinetic energy occurs near the bottom of the flux tube and it steadily decreases with 
altitude; near the apex the ions are super cold with an average energy less than that of the ambient 
0 + ions in the ionosphere. Such cold ions may contribute to the hidden cold core, which are 
difficult to measure using particle detectors, but affect the characteristic plasma frequency. Such 
cold thermal 0 + ions have been observed up to altitudes of several thousand kilometers, and so far 
no explanation has been found as to how they reach such heights [Abe et al, 1993; Peterson et al, 
1993]. Our study reveals that a low level of transverse ion heating at low altitudes in conjunction 
with the gravity creates a cold 0 + population extending to high altitudes. 

The above features of 0 + produced by transverse heating occurs only for sufficiently low 
heating. For strong heating producing 0 + ions greatly in excess of 10 eV energy, most ions 
escape and the features discussed above become quite weak. 
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3.3 0 + Flux 

Outilux: In the past several studies have been performed on 0 + outflux [ Barakat and Schunk, 
1983; Li et al, 1988; Korosmezey et al , 1992; Wilson , 1993], We show here that transverse ion 
heating lasting over a few minutes produces a time varying flux; the evolution of the flux as a 
function of altitude and time is examined here. For the low heating rate discussed in this paper, 
the evolution of outflux of 0 + at early times is shown in Figure 11a. In this figure we have 
plotted normalized outflux r = n 0+ u 0 r 3 /(R^, + 1000 km) 2 \ the normalization factor 
r o = n a Vfl , which is 5 • 5 x 10 9 ions cm~ 2 s~ x , if we assume n a = 10 4 cm -3 and T 0 = 0.3 eV . 
We find that when the heating occurs at geometric distance r >S3\0 km ( h> 1940 km), a 

large flux develops even below the heating region and it extends all the way down to the base at 
r =r b , which has the height of 980 km. One minute after the heating begins, the maximum 
outflux at the base is r = 6x\0~ 2 r o . If we assume n a = 10 4 cm~ 3 , and T o =0.3eV, 
r o - 5- 5 x 10 9 ions cm~ 2 s~ l and r = 3 x 10 8 ions cm~ 2 s ~ l . This large flux develops due to 
average upward flow developing all the way down to the base. Figure 1 lb shows the evolution of 
the average flow velocity giving rise to the flux. At / = 0 the flow velocity is zero or less than 
zero all along the flux tube. At t = 1 min., the flow velocity at the base is 0.06 s 324 mis, 

and decreases with the increasing altitude. At later times, when the heated ions move upward, the 
flow velocity increases at higher altitudes while it decreases at lower altitudes and eventually it 
approaches zero there. This is illustrated by the curve for t = 6 min. in Figure 1 lb. It is 

interesting to examine how such a large flux develops at altitudes below the heating region. It is 
not due to any upward acceleration of ions below the heating altitude. Before the heating, 
downward and upward fluxes are nearly balanced giving a zero flux, but when heating occurs 
some downflowing ions in the initial population are not able to flow downward any more because 
of the mirror force. This creates an imbalance in favor of an upward flux. 

Figures 12a and 12b show the further evolution of the outflux by plotting F at some 
selected times, at t = 10, 30 and 60 mins, in Fig. 12a and at t = 2,3 and 4 hours in Fig. 12b. We 

see that as the outflux spreads upward its maximum value decreases; at 


15 



t = 30 min., /' max = 2-25x10 J = 1 • 2 x 1 0 ' ions cm ~s ! ana for t> 2 hrs. the flux is less 

*3 r ^ i 

than 10 i 0 =5-5x10 cm ~s . The above values of the outflux should be compared against 
the flux produced by warm electrons, when electron temperature is enhanced from 3000° k to 
10,000% the outflux increases from a negligibly smail value to 0-4x10 1 ions cm' 2 / s~ ] 
[Barakat and Schunk, 1983], A similar magnitude of flux in response to elevated 0 + temperature 
was reported by Li et al [1988], We see that even a weak transverse ion heating can produce a 
substantially large flux extending to r ~3x\Q 4 km in comparison to that produced by an 

enhancement in the electron and ion temperatures. However, the transverse heating, occurring 
over a brief time period of a few minutes, produces oniv a transient outflux. Figure 12b also 
shows that relatively weak outflux developing at reiativeiv iow altitudes at later times. For 
example, note the outflux below r - 15,000 km at all three times shown in the figure. Such fluxes 

are the consequences of secondary expansion fronts which develop when the down-flowing ions 
mirror upward. The pulse of flux just below r = 1 5,000 km at / = 240 min. corresponds to the 
expansion front shown by the arrow in Figure 4a. The three peak outfluxes for r < 1 5,000 km are 
approximately 5-5 xlO 6 , 1-lxlO 6 , and 5 • 5 x 1 0 5 ions cm ~ 2 s~ 1 at t = 2, 3 and 4 hrs., 
respectively. 

Downtlux: From our previous discussion on down moving ions, we should expect downward 
fluxes developing in response to transverse ion heating. We discuss here the magnitude of the 
downward flux. We find two separate instances of downward flux. The first case occurs at very 
early stage of the heating and confined to quite low altitudes, even below the altitudes where 
heating occurs. The second case occurs after the upward moving ions are able to travel a 
substantial distance and then they are pulled down by the gravity. 

Figure 13 shows the low-altitude downward flux occurring at early times. At / = 10 min., 
the maximum downward flux is 8x10 r o = 4-4x10 ions cm ~s right near the base of the 
flux tube. How is this large downward flux created? In order to answer this question, we note 
from Figure 11a that a large upward flux develops at eariv time even at altitudes below the 
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heating. As noted earlier, this large outflux is due to the imbalance between the downward and 
upward flowing ions below the heating region caused by the heating of some of the downflowing 
ions at higher altitudes. This imbalance favoring upward flux does not last long. Since ions 
carrying the additional upward flux are not energized, they do not travel far and are turned 
downward by the gravity, creating a substantial downward flux at low altitude at an early time; for 
example, at / = 10 min. the downward flux is limited to altitude below the heating region for 
r > 8310 km . At later time the downward flux is seen to penetrate into the heating region. This 
depends on the initial energies of the up-flowing ions carrying the upward unbalanced flux and the 
perpendicular energy they may receive in the heating region. If they receive relatively large 
perpendicular energy, they may not return and this reduces the downward flux as seen for t - 30 
and 60 mins, in Fig. 13. 

The large initial downward flux at low altitudes gives way to smaller downward flux when 
ions traveling to higher altitudes begin to return. Figure 14a shows such downward fluxes for 
t - 2, 3, and 4 hrs. In order to understand the downward flux, we have plotted the average flow 

velocity of 0 + in Figure 14b at the same times as in Figure 14a. At / = 2 hrs., we expect 
downward flux over the altitude range 10,000 km <r < 19,000 km , for which the flow velocity 

V < 0. We see the downward flux for this time in Figure 14a; the maximum flux occurs at about 
r = 1 1,000 km and its magnitude is 1-6 x 1 0 6 ions cm~ 2 s~ x . At later times the downward flux 

generally decreases, but it extends to higher altitudes, as we discussed earlier in connection with 
the phase space plots in Figures 3 and 4. We notice from Figure 14a that there are dropouts in 
the downward flux in some altitude intervals. Such dropouts occur at altitudes, where upward 
expanding density fronts are located at a given time. In the expansion front the flux is upward as 
shown in Figure 12b. 

In view of Figure 14b, we may further amplify our discussion on the drift velocity of the 
downflowing ions. We notice that at times later than 2 hrs., after which expansion fronts begin to 
appear, the maximum downward drift velocity occurs just above the expansion front and the 
height where it occurs increases with increasing time. Furthermore, its magnitude increases with 
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the altitude. This increase in the magnitude of the maximum drift velocity is simply the 
consequence of the relatively more energetic ions traveling to higher heights and in the process of 
falling require a higher downward drift velocity. 

4. Conclusion and Discussion 

In this paper we have studied the time response of 0 + flow to a low level of transverse 
heating in the polar ionosphere. The heating is allowed to occur over a few minutes and the 
heating rate is kept sufficiently low so that the bulk of the heated ions have a few eV energies. 
Such ions do not escape the downward gravitational pull. Therefore, the heated ions first move 
upward as they are acted upon by the mirror force, but they are eventually turned downward by 
the gravity. While moving downward, the ions mirror back eventually and execute a bounce 
motion in the flux tube. The dynamics and the fluxes of 0 + associated with the initial upflow in 
response to the transverse heating and with the subsequent bounce motion are examined in this 
paper. The main results of this paper are as follows: 

1) When a low level transverse ion heating occurs in the topside ionosphere, it generates a large 
outflow of 0 + even at altitudes below the heating region. This outflux is the consequence of 
the imbalance between the upward and downward fluxes. A large outflux below the heating 
region lasts over a few minutes, after which a large downward flux develops lasting again just 
over a few minutes. The interesting fact is that ions participating in such an event of large 
upward flux followed by a large downward flux below the heating region are generally 
unenergized. 

2) The initial upward expansion of the heated ions significantly enhances 0 + density to 
geocentric distance up to several Re . However, the ions transported to high altitudes do not 
stay there. They are pulled down by the gravity. The height reached depends on the energy 
gained during the heating period. The lower the energy gain, the lower is the height reached 
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by the ions. Thus, during the course of the evolution of the heated ions it is possible to have 
up-flowing ions at high altitudes and down-flowing ions below them. 

3) The average flow velocity of the down-flowing ions maximizes above the topside ionosphere 
because the dense cold plasma inside the ionosphere reduces the average velocity of the ions 
accelerating downward under the influence of the gravity even before the mirror force begins 
to dominate the gravity. The downward flow velocity of 0 + seen from the simulation is 
comparable to those observed from DE 1 [ Chandler , 1995]. 

4) When downward flowing ions mirror back, they generate a series of expansion fronts moving 
upward from the ionosphere. The formation of such expansion fronts is seen to last over 
several hours in response to a brief heating time of just a few minutes. The width of the 
expansion front is determined by the heating time and the spread in the 0 + parallel velocity, 
which evolves due to the action of the mirror force and the gravity. Figure 6 shows that the 
expansion front is relatively thick at lower altitudes and becomes thinner as it moves upward 
because the 0 + ions progressively cool. 

5) In response to the formation, and upward propagation of the expansion fronts, the local 0 + 
densities show nearly periodic oscillations with a time period increasing with altitude. The 
period ranges from about an hour at the topside ionosphere to several hours at the altitude of 
a few Re. 

6) The formation of expansion fronts ceases when the heated ions, having a range of initial 
energies, disperse in space. This happens in about one day. 

7) When the trapped heated ions have dispersed in space, the resulting 0 + density profile extends 
to large distances with significantly enhanced 0 + densities. Thus a low level transverse 
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heating occurring in the topside ionosphere is a possible mechanism for the supply of 0 + ion 
to the distant part of the magnetosphere. 

8) When the trapped ions have dispersed in space, the ions with largest kinetic energies occur in 
the topside ionosphere, while the coldest ions with energies less than the characteristic 
ionospheric energy ( kT a ) occur at large geocentric distances near the apex of the trapped ion 

trajectories. Thus a low-altitude heating is a mechanism for producing the cold plasma 
population. Such populations are difficult to measure by particle detectors while they affect 
the plasma frequency. This yields a higher estimate of plasma density from wave 
measurements compared to that from direct particle measurements. However, cold thermal 
ions have been observed at distances of several tens of thousand of kilometers from DE-1 
[Abe et al, 1993; Peterson , 1993], but without a plausible explanation. Our study here 
provides a plausible explanation in terms of a low level transverse heating in the topside 
ionosphere in conjunction with the downward gravitational force, which reduces the kinetic 
energy as the ions travel to higher altitudes. 

9) Simulation predicts a range of fluxes for both up- and down-flowing ions. The 0 + outfluxes 
resulting from the transverse heating are generally transient. The maximum magnitude of such 
outfluxes is about ~ 10 ions cm ~s , as predicted by previous steady-state models [Barakat 

and Schunk, 1983; Li et al, 1988], In addition to the upward flux of the heated ions, a 
transient heating event causes a relatively large transient outflux of 0 + ions below the heating 
region. This is simply caused by the imbalance in the upward and downward fluxes when the 
downward flowing ions mirror back due to the transverse heating. 

10) Except for a large downward flux during the early stage of the heating at low altitudes below 
the heating region, the downward flux of 0 + above the dense 0 + plasma is found to be 
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~ 10 ions cm s . The downward drift and flux show a great deal of variability with both 
time and altitude. 

Finally, we point out that in the one-dimensional calculations presented here, the 
convection electric field was ignored. Therefore, the resuits presented here are the response of a 
single flux tube as it convects, keeping all its plasma contents as under the frozen-in-field 
approximation. In other words, an observer moving with the convecting flux tube should see the 
processes described here. In the rest frame of reference, the inclusion of the convection electric 
field gives rise to a hopping motion as revealed by a two-dimensional trajectory calculation of 
Horwitz [1984], In order to assess the influence of the hopping motion on the dynamical behavior 
of the superthermal 0 + ions, a two-dimensional model based on self-consistent electric field is 
needed. However, for a heating region extending both in longitude and latitude, the results 
presented here should be qualitatively applicable; in such a situation the rising and falling ions 
should uniformly repeat themselves. The most important result of this paper on the creation of a 
cold 0 + population at high altitudes should not be modified by the convection electric field. A 
multitude of heating events randomly distributed in space and time should produce the cold 0 + 
ion population at altitudes up to several thousand kilometers in a large volume of space covering 
the polar and auroral regions. 
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FIGURE CAPTIONS 


Fig. 1. Geometry of the flux tube. We consider a flux tube from a geocentric distance r b to 
'max ■ 1° the simulation macroparticles representing 0 + and H + ions are injected into the flux tube 
at the base r - r b and those exiting the flux tube at r = are assumed to be lost. 

Fig. 2. The properties of the ambient polar wind (t = 0) . (a)-(b) Phase space plots of 0 + ions. 
(c)-(d) Phase space plots for H + ions. Left and right hand panels for the phase space plots are in 
r-V\ | and r-V L planes, respectively. Each dot in the phase space represents a macro 0 + or H + 
ion. (e) 0 + and H + density, and electric field profiles, (f) 0 + and H + average flow velocities 
and H + outflux F . 

Fig. 3. Same as Figure 2, but for later times, (a)-(b) t = 5 min., (c)-(d) t - 20 min., (e)-(f) 
f = 60 min. , and (g)-(h) / = 120 min. Note that at t = 120 min. , both down and upflowing ions 
occur. 

Fig. 4. Same as Figure 3, but for later times, (a) -(b) t - 4 hrs., (c) -(d) t = 6 hrs., (e) - (f) 
t = 8 hrs., and (g) -(h) / = 12 hrs. Note the rising 0 + clouds as indicated by arrows. 

Fig. 5. Density profiles during the initial upward expansion of the heated ions. The solid-line 
curve shows the 0 + density profile for the ambient polar wind ( t = 0 ). Note the beginning of re- 
expansion at / = 180 min. The inset shows the temporal evolution of the electric field at 
r = 9230 km . 

Fig. 6. Same as Fig. 5, but for later times showing further motion of the expansion front 
appearing at / = 180 min. Note a new expansion front appearing at / = 360 min. 
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Fig. 7. Average flow velocity associated with the expansion front shown in Figure 6. 

Fig. 8. Further upward motion of the new expansion front at t = 360 min. , as noted in Fig. 6. 
The new expansion front, indicated by short arrows, has reached r = 16,000 km at / = 8 hrs. 
The old expansion front (indicated by long arrows) has reached r 5 23,00 km at the same time. 

Fig. 9. Characteristic density oscillations at some selected geocentric distances in response to the 
formation and propagation of fronts shown in Figs. 6 to 8 

Fig. 10. Density profile associated with the trapped superthermal ions at / = 24 hrs. The initial 
density profile without the heating is also shown. Note the enhancement in the density for 
r > 10,000 km . The inset shows the phase-space plots for the diffused superthermal ions trapped 
inside the flux tube. Note that the wannest population of 0 + occurs near the bottom of the flux 
tube, while near the apex of the trapped 0 + ions are quite cold. 

Fig. 11. (a) Outflux at early times, and (b) average flow velocity associated with the outflux in 

(a)- 

Fig. 12. Outflux at times later than were shown in Figure 1 la: (a) Fluxes at t - 10, 30 and 60 
mins., and (b) Fluxes at / = 2, 3 and 4 hrs. Note the "pulse" of outflux developing for 
r < 15,000 km ; the flux is generated by the expansion fronts shown in Figure 6. 

Fig. 13. Downward flux at relatively early times. 

Fig. 14. (a) Downward flux at times later than those shown in Figure 13, and (b) average flow 
velocity associated with the downward flux in (a). 
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Abstract: 

Normally the gravitationally bound heavy 0 + ions in the Earth's ionosphere are in a 
diffusive equilibrium. However, when energized to superthermai energies of a few eV transverse 
to the geomagnetic field, the combined effects of the downward gravitational and the upward 
electric and mirror forces produce interesting flow patterns in the vertical direction like in a 
pulsating fountain. This flow pattern is studied by means of a particle-in-cell code. 

It is well known that the topside ionosphere is commonly populated with a variety of ion 
species such as H + , 0 + , and H* . The light ions (H + and H*) readily escape the downward 
gravitational force even when they have a small characteristic energy of about a few tenths of an 
eV. On the other hand, for the heavy 0 + ions the escape requires some sort of energization. For 
example, an 0 + ion escapes when its energy W>W es =m Q +g 0 R%!r, where m QJr is 0 + 
mass, g a is the Earth's gravity at its surface, R e is the Earth's radius, and r is the geocentric 
distance in the topside ionosphere. Using R e =6371 km, g Q = 9.8 ms~ 2 and r = 7350 km, 
the escape energy W ea = 7.5 eV . For W < W es the 0 + ions can be transported to large altitudes, 
but they eventually fall down. If the energization occurs in the degree of freedom parallel to the 
magnetic field, the falling ions are lost into the ionosphere. On the other hand, when the heating 
occurs in the degree of freedom transverse to the geomagnetic field, the 0 + ions are trapped in 
the flux tube and undergo a bounce motion under the influence of the downward gravitational 
force and the upward mirror force associated with the downward gradient in the geomagnetic 
field B . The upward polarization electric field associated with the density gradient also plays a 
critical role in the bounce motion. 




We study vertical plasma flow in geomagnetic open flux tubes, which are typical of the 
polar ionosphere. The topside ionospheric base is assumed to be at r -r b =7350 km. We 
simulate the flux tube from r -r b to a geocentric distance r - r mnv = 31,350 km. In order to 
study the plasma flow, we employ a particle code [1], We inject macroparticles corresponding to 
0 + and H 4 " ions into the flux tube at the boundary at r = r b and subsequently follow their 
dynamics under the influence of the gravity, polarization electric field and the mirror force 
associated with the gradient in the geomagnetic field. The injected particles are chosen from 
Maxwellian velocity distributions with a zero drift for 0 + and a drift of 2V A for H 4 ", where 
is the thermal velocity corresponding to the temperature T a of H 4- ions at the base. We assume 
that O 4 " and electron temperatures at r = r b are also T 0 . When a steady-state polar wind flow is 
set up in the flux tube, we initiate the transverse heating event for O 4 * ions at geocentric distances 
r > 83 1 1 km . The heating at such altitudes is expected to be caused by wave-particle interactions 
and it is incorporated in a random fashion as described by Brown et al [2]. We present here the 
flow pattern of 0 + ions. The H 4 " ions do not show any interesting feature except for their 
supersonic escape. Since H* ions have extremely low concentration compared to that of O 4 * and 
H 4 ", they are not included in the model. 

We show the bounce motion of the heated ions and its temporal evolution through a 
sequence of phase space plots after a brief heating event of 5 minutes at an average rate of about 
0.04 k T 0 Is. When kT 0 =Q3eV, this leads to a maximum energization of about 15 eV. 
Therefore, some of the ions, having energies W > W u escape while those with W < W es are 
trapped. Figs. 2a, 2b and 2c show the initial diffusive state of the 0 + . Note that the left, middle 
and right panels in Fig. 1 are the phase space plots in r - V\\ , r - and r-W planes, where r 
is the geocentric vertical distance, and V ± are the 0 + velocity components parallel and 
perpendicular to the geomagnetic field line, respectively, and IF = }m 0 (Fj, 2 +{^), that is, it is 
the total kinetic energy. Each dot in the figures represent an 0 + ion. Figures 2a, 2b and 2c show 
that the bulk of O 4 " ions are generally confined below an altitude r = 10 4 km . 


The state of the 0 + flow after 6 hours of the brief heating period is shown in Figs 2d, 2e 
and 2f. Note that some ions have reached near the top of the simulated flux tube and have 
V\\ > 0 , implying that they will escape. While those having V\\ < 0 are falling. Those ions are 
transported to high altitudes by the mirror force F m x - L m Q V\B~^ ( dB / d r). The transverse 
heating increases Fj_ and enables the upward flow. But the gravity pulls them downward if they 
do not have sufficiently large energy to escape. 

The fallin g ions mirror back and expand upward and they appear in clumps; the clumps are 
seen in the form of evolving bands in the r-V\ \ plane as seen from Figs 2d, 2g, and 2j. The 
bands give the appearance of a water fountain. The bands can be seen as clusters of ions in the 
r-V j_ plane. V L decreases with increasing r due to the conservation of the magnetic moment 
M =2 IB ■ 1° terms °f energy, the ion clusters appear in the form of "bananas," which 

simply manifest that the upgoing ions have slightly more kinetic energy than the downflowing 
ions. In both r~V ± and r-W planes, the upflowing > 0) ions are shown in green while 
the downflowing are shown in blue. 

The cluster formation occurs by the modulation in the electric field E due to the falling 
and rising 0 + ions because E depends on the density gradients. This modulation of E also 
causes the difference in the kinetic energies of the upgoing and downgoing ions in any given 
cluster. 
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Figure Captions 


Fig. 1 . Phase space plots in r- Vjj (left panels), r-V ± (middle panels), and r-W (right 
panels) at several times are shown. Note that the velocity components are normalized with 
respect to the thermal velocity of H + ions which is about 5.4 km! s for kT a - 0.3 eV and 
energy is normalized in terms of k T g . In the middle and right panels, blue represents falling 0 + 
ions, while green represents upflowing ones. 



